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I Abstract 

'^ We consider the low energy implications including particle spectroscopy of SO (10) 

inspired t-b-r Yukawa coupling unification with /x < 0, where /x is the coefficient of the 
^ bilinear Higgs mixing term of the minimal supersymmetric standard model (MSSM). 

We employ non-universal MSSM gaugino masses induced by SO(IO) invariant dimen- 
sion five operators, such that the total number of fundamental parameters is precisely 
the same as in Yukawa unified supersymmetric SO(IO) models with universal gaugino 
Ph masses and n > 0. We find that t-b-r Yukawa unification with yU < is compatible 

p with the current experimental bounds, including the WMAP bound on neutralino 

^ dark matter and the measured value of the muon anomalous magnetic moment. We 

I— I present a variety of benchmark points which include relatively light squarks (~ TeV) 

I of the first two families and an example in which the bottom and top squarks are 

>• lighter than the gluino. This is quite distinct from Yukawa unification with /i > 

^R and universal gaugino masses in which the gluino is the lightest colored sparticle and 

psj the sqaurks of the first two families have masses in the multi-TeV range. 
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1 Introduction 

Supersymmetric (SUSY) 5*0(10) grand unified theory (GUT), in contrast to its non- 
SUSY version, yields third family (t-b-r) Yukawa unification via the unique renor- 
malizable Yukawa coupling 16 ■ 16 ■ 10, if the Higgs 10-plet is assumed to contain 
the two Higgs doublets Hu and H^ of the minimal supersymmetric standard model 
(MSSM) [1]. The matter 16-plet contains the 15 chiral superfields of MSSM as well 
as the right handed neutrino superfield. The implications of this Yukawa unification 
condition at Mq ~ 2 x lO'^^ GeV have been extensively explored over the years [1, 2]. 
In S'O(IO) Yukawa unification with /i > and universal gaugino masses, the gluino 
is the lightest colored sparticle [3, 4], which will be tested [5] at the Large Hadron 
Collider (LHC). The squarks and sleptons, especially those from the first two families, 
turn out to have masses in the multi-TeV range. Moreover, it is argued in [3, 4] that 
the lightest neutralino is not a viable cold dark matter candidate in 5*0(10) Yukawa 
unification with /i > and universal gaugino masses at Mq 

Spurred by these developments we have investigated t-b-T Yukawa unification [4, 
6, 7] in the framework of SUSY 5f/(4), x SU{2)l x SU{2)r [8] (4-2-2, for short). The 
4-2-2 structure allows us to consider non-universal gaugino masses while preserving 
Yukawa unification. An important conclusion reached in [4, 6] is that with same 
sign non-universal gaugino soft terms, Yukawa unification in 4-2-2 is compatible with 
neutralino dark matter, with gluino co-annihilation [4, 5, 6, 9] being a unique dark 
matter scenario for /i > 0. 

By considering opposite sign gauginos with /i < 0, M2 < 0, M3 > (where fi is 
the coefficient of the bilinear Higgs mixing term, M2 and M3 are the soft supersym- 
metry breaking (SSB) gaugino mass terms corresponding respectively to SU{2)l and 
S'f7(3)c). It is shown in in [7] that Yukawa coupling unification consistent with the 
experimental constraints can be implemented in 4-2-2. With /i < and opposite sign 
gauginos, Yukawa coupling unification is achieved for mo > 300 GeV, as opposed to 
niQ > 8 TeV for the case of same sign gauginos. The finite corrections to the b-quark 
mass play an important role here [7]. By considering gauginos with M2 < 0, M3 > 
and /i < 0, we can obtain the correct sign for the desired contribution to {g — 2)^ [10]. 
This enables us to simultaneously satisfy the requirements of t-b-r Yukawa unification 
in 4-2-2, neutralino dark matter and {g — 2)^, as well as a variety of other bounds. 

Encouraged by the abundance of solutions and coannihilation channels available 
in the case of Yukawa unified 4-2-2 with M2 < and yU < 0, it seems natural to 
explore Yukawa unification in SO(IO) GUT (with M2 < and /x < 0). It has been 
pointed out [11] that non-universal MSSM gaugino masses at Mq can arise from 
non-singlet F-terms, compatible with the underlying GUT symmetry such as SU(5) 
and SO(IO). The SSB gaugino masses in supergravity [12] can arise, say, from the 



following dimension five operator: 

T7ab 



2Mi 



PI 



A"A'' + c.c. (1) 



Here A° is the two-component gaugino field, F"''' denotes the F-component of the 
field which breaks SUSY, the indices a, b run over the adjoint representation of the 
gauge group, and Mpi = 2.4 x 10^^ GeV is the reduced Planck mass. The resulting 
gaugino mass matrix is {F°'^)/M-p\ where the supersymmetry breaking parameter 
{F"-^) transforms as a singlet under the MSSM gauge group SU{?>)c x SU{2)l x U{1)y- 
The F"-^ fields belong to an irreducible representation in the symmetric part of the 
direct product of the adjoint representation of the unified group. 
In SO(IO), for example, 

(45 X 45)5 = 1 + 54 + 210 + 770 (2) 

If F transforms as a 54 or 210 dimensional representation of SO(IO) [11], one obtains 
the following relation among the MSSM gaugino masses at Mq '■ 

M3 : M2 : Ml = 2 : -3 : -1, (3) 

where Mi,M2,M3 denote the gaugino masses of f/(l), SU{2)l and S'f/(3)c respec- 
tively. The low energy implications of this relation have recently been investigated 
in [13] without imposing Yukawa unification. 

The outline for the rest of the paper is as follows. In Section 2 we summarize 
the scanning procedure and the experimental constraints that we have employed. In 
Section 3 we present the results from our scan and highlight some of the predictions 
of an SO(IO) model with /i < and the non-universal MSSM gaugino masses at Mq 
related by Eq.(3). We display some benchmark points which can be tested at the 
LHC. Our conclusions are summarized in Section 4. 



2 Phenomenological Constraints and Scanning Pro- 
cedure 

We employ the ISAJET 7.80 package [14] to perform random scans over the fun- 
damental parameter space. In this package, the weak scale values of gauge and 
third generation Yukawa couplings are evolved to Mq via the MSSM renormalization 
group equations (RGEs) in the DR regularization scheme. We do not strictly enforce 
the unification condition gs = gi = (72 at Mq, since a few percent deviation from 
unification can be assigned to unknown GUT-scale threshold corrections [15]. The 
deviation between gi = g2 and g^ at Mq is no worse than 3 — 4%. For simplicity we 



do not include the Dirac neutrino Yukawa coupling in the RGEs, whose contribution 
is expected to be small. 

The various boundary conditions are imposed at Mq and all the SSB parameters, 
along with the gauge and Yukawa couplings, are evolved back to the weak scale 
Ml. In the evaluation of Yukawa couplings the SUSY threshold corrections [16] are 
taken into account at the common scale Mgusy = y/^^W^'^H^- The entire parameter 
set is iteratively run between Mz and Mq using the full 2-loop RGEs until a stable 
solution is obtained. To better account for leading-log corrections, one-loop step-beta 
functions are adopted for gauge and Yukawa couplings, and the SSB parameters irti 
are extracted from RGEs at multiple scales mi = ini{ini). The RGE-improved 1-loop 
effective potential is minimized at MgusY? which effectively accounts for the leading 
2-loop corrections. Full 1-loop radiative corrections are incorporated for all sparticle 
masses. 

The requirement of radiative electroweak symmetry breaking (REWSB) imposes 
an important theoretical constraint on the parameter space. In order to reconcile 
REWSB with Yukawa unification, the MSSM Higgs soft supersymmetry breaking 
(SSB) masses should be split in such way that m^ /""^h^ > 1-2 at Mq [17]. As 
mentioned above, the MSSM doublets reside in the 10 dimensional representation 
of SO(IO) GUT for Yukawa unification condition to hold. In the gravity mediated 
supersymmetry breaking scenario [12] the required splitting in the Higgs sector can 
be generated by involving additional Higgs fields [18], or via D-term contributions 
[19]. Another important constraint comes from limits on the cosmological abundance 
of stable charged particles [20] . This excludes regions in the parameter space where 
charged SUSY particles, such as fi or ti, become the lightest supersymmetric particle 
(LSP). We accept only those solutions for which one of the neutralinos is the LSP 
and saturates the WMAP bound on relic dark matter abundance. 

We have performed random scans for the following parameter range: 

< mo < 5 TeV 

< niH^ < 5 TeV 

< mu^ < 5 TeV 

< Ml/2 < 2 TeV 

35 < tan/3 < 55 

-3 < Aq/tjio < 3 (4) 

with /i < and mt = 173.1 GeV [21]. Note that our results are not too sensitive to 
one or two sigma variation in the value of mt [18]. We use mb{m,z) = 2.83 GeV which 
is hard-coded into ISAJET. The set of parameters presented above is usually referred 
to as NUHM2 [22]. This choice of parameter space was informed by our previous 
works on t-b-r Yukawa Unification [6, 18]. 



Employing the boundary condition from Eq.(3) one can define the MSSM gaugino 
masses at Mq in terms of the mass parameter M1/2 : 



Ml = 


-Ml/2 


M2 = - 


-3Mi/2 


M3 = 


= 2M1/2 



(5) 

In scanning the parameter space, we employ the Metropolis-Hastings algorithm 
as described in [23]. The data points collected all satisfy the requirement of REWSB, 
with the neutralino in each case being the LSP. After collecting the data, we impose 
the mass bounds on all the particles [20] and use the IsaTools package [24] to imple- 
ment the various phenomenological constraints. We successively apply the following 
experimental constraints on the data that we acquire from ISAJET: 

ruh (lightest Higgs mass) > 114.4 GeV [25] 

BR{B, -^ H+H-) < 5.8 X 10-s [26] 

2.85 X 10-^ < BR{b -^ 57) < 4.24 x 10"^ {2a) [27] 

^■^^ < ""mCZiTsT < 2.41 (3a) [27] 

QcDuh'' = 0.111+°:0f7 (5a) [28] 

0<A(^- 2)^/2 < 55.6x10-1° [10] 



3 Yukawa Unification and Particle Spectroscopy 

We next present the results of the scan over the parameter space listed in Eq.(4). 
In Fig. 1 we show the results in the R — rriQ, R — tan/3, R — M1/2 and M1/2 — ""^o 
planes. Gray points are consistent with REWSB and neutralino LSP. Green points 
satisfy particle mass bounds and constraints from BR{Bs -^ fi^^~), BR{h -^ 57) 
and BR{Bu — )■ rur). In addition, we require that green points do no worse than the 
SM in terms of {g — 2)^. Orange points belong to a subset of green points and satisfy 
the WMAP bounds on Xi dark matter abundance. In the M1/2 — '"^o plane, points in 
brown represent a subset of yellow points that are consistent with Yukawa coupling 
unification to within 10%. 

In the R - rriQ plane of Fig. 1 we see that with both /i < and M2 < 0, we 
can realize Yukawa unification consistent with all constraints mentioned in Section 2 
including the one from [g — 2)^. This is possible because for /x < 0, we can implement 
Yukawa unification for relatively small mo(~ 500 GeV), and, in turn, {g — 2)^ obtains 
the desired SUSY contribution proportional to /XM2. This is more than an order of 
magnitude reduction on the mo values required for Yukawa unification with /i > and 
universal gaugino masses. In the present with 10% or better t-h-r Yukawa unification 
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Figure 1: Plots in R — rrio, R — tan/3, R — M1/2 and M1/2 — mo planes. Gray points 
are consistent with REWSB and neutralino LSP. Green points satisfy particle mass 
bounds and constraints from BR{Bs — )■ fi^fi^), BR{h — )• 57) and BR{Bu — )• tz/,-). In 
addition, we require that green points do no worse than the SM in terms of {g — 2)^. 
Orange points belong to a subset of green points and satisfy the WMAP bounds on 
x\ dark matter abundance. In the M1/2 — mo plane, points in brown represent a 
subset of yellow points and satisfy Yukawa coupling unification to within 10%. 



we obtain a relaxation also of Aq values similar to the 5'f/(4)c x SU{2)l x SU{2)ji 
model in [7], with —2.5 < Ao/mo < 2. Our observation about relaxing the possible 
range of tan (3 that accommodates Yukawa unified models is explicitly shown in the 
R - tan /3 plane. In the R - mi/2 plane of Fig. 1 we see that employing the boundary 
conditions for gauginos presented in Eq. (5), the lightest neutrahno mass can be as 
low as 15 GeV consistent with all constraints mentioned in Section 2 including the 
one from {g — 2)^. Note that it is impossible to realize a neutralino mass as this in the 
universal gaugino case due to the chargino mass constraint. A narrow orange strip 
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Figure 2: Plots in the rrif - m^o, itla - m^o, 
The gray points satisfy the requirements of REWSB and Xi LSP. Dark green points 
satisfy particle mass bounds and constraints from BR{Bs — ?■ fi^fi"), BR{b — )■ S7) 
and BR{Bu — ?■ rur). In addition, we require that green points do no worse than the 
SM in terms of {g — 2)^. Light green points are a subset of these points which also 
satisfy Yukawa unification. We show in the ruf - m^o and m-± - m^o planes the unit 

•' ' Al Xl Xl -^ 

slope lines representing the respective coannihilation channels. In the m^ - m^o plane 
we show the line niA = 2Tn^o that signifies the A resonance channel. 



for low Ml/2 values indicates the existence of Z and light Higgs resonance solutions 
for neutralino dark matter. Actually, there are two very narrow strips, one around 
45 GeV and a second around 60 GeV, even though they appear as one strip in the 
figure. In order to better visualize the magnitude of the sparticle masses consistent 
with t-b-T Yukawa unification, we present our results in the M1/2 — m-o plane, where 
the brown points correspond to Yukawa unification better than 10%. 

In Fig. 2 we show the relic density channels consistent with Yukawa unification 
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Table 1: Sparticle and Higgs masses (in GeV). All of these benchmark points sat- 
isfy the various constraints mentioned in Section 2 and are compatible with Yukawa 
unification. Point 1 depicts a solution corresponding to the A funnel region. Points 
2 and 3 display the light Higgs and Z-resonance solutions, while point 4 represents 
the stau coannihilation solution. 



in the m^ - m^o, itla - m^o, m-± - fn^o, and rriq - rrig planes. Gray points shown in 
this figure satisfy the requirements of REWSB and y^ LSP. Dark green points satisfy 
the particle mass bounds and constraints from BR{Bs — ?■ fi^fi^), BR{b — )■ S7) and 
BR{Bu — )■ TUr). the green points do no worse than the SM in terms of {g — 2)^. 
The light green points represent a subset of the dark green points, and correspond 
to 10% or better t-b-T Yukawa unification. This choice of color coding is infiuenced 
from displaying the sparticle spectrum with and without t-b-r Yukawa unification, 
while still focussing on all the other experimental constraints. The idea is to show the 
myriad of solutions that implement Yukawa unification and are consistent with all 
known experimental bounds except that on relic dark matter density from WMAP. 
The appearance of a variety of Yukawa unified solutions with a very rich sparticle 
spectrum is a characteristic feature of /x < [7] . 

We can see in Fig. 2 that a variety of coannihilation and resonance scenarios are 
compatible with Yukawa unification and neutralino dark matter. Included in the rriA 
- m^o plane is the line niA = 2m^o which shows that the A-funnel region is compatible 
with Yukawa unification. In the ruf - m^o plane in Fig. 2, we draw the unit slope 
line which indicates the presence of stau coannihilation scenarios. From the m-± - 

Xl 

m^o plane, it is easy to recognize the light Higgs (h) and Z resonance channels. We 
expect that other coannihilation channels like the stop coannihilation scenario are 
also consistent with Yukawa unification, although we have not found them, perhaps 
due to lack of statistics. 

Let us remark on the low mass neutralino solutions that we have found in our 
model (with /x < 0). Because of the Mq scale gaugino mass relations in Eq.(5), it is 
possible in principle to have small Mi values, thus giving rise to a light neutralino. 
The neutralino mass nonetheless is bounded from below because of the relic density 
bounds on dark matter. The SO(IO) model with non-universal gaugino masses, as in 
this paper, has all the ingredients to bring down the neutralino mass to the lowest 
possible value consistent with the various constraints. The solution with the neu- 
tralino (mass ~ 43 GeV) is consistent with Yukawa unification and corresponds the 
Z-resonance dark matter scenario. 

Finally, in Table 1 we present some benchmark points for the SO(IO) t-b-T Yukawa 
unified model with /i < and non-universal gaugino masses. All of these points 
contain WMAP compatible with neutralino dark matter and satisfy the constraints 
mentioned in Section 2. Point 1 depicts a solution with essentially perfect Yukawa 
unification corresponding to the A funnel region. Points 2 and 3 correspond to the 
light Higgs and Z-resonance solutions, while point 4 represents the stau coannihilation 
solution. It is interesting to note that for the light Higgs and Z-resonance solutions, 
there is an upper bound on the gluino mass. Employing the boundary condition in 
Eq.(5) this turns out to be nig ^ 700 GeV. Hence, the light Higgs and Z-resonance 
solutions are not compatible with this model if the gluinos are founded to be heavier 



than ^ 700 GeV. 

4 Conclusion 

We have shown that SO(IO) t-b-r Yukawa unification with /i < and non-universal 
gaugino masses is nicely consistent with all available experimental data. We have 
considered a variety of WMAP compatible neutralino dark matter scenarios, including 
some examples in which the LSP neutralino can be rather light, about half the Z- 
boson mass (Z-resonance solution) or the SM-like Higgs mass (light Higgs resonance 
solution). Neutralino dark matter solutions corresponding to the A- funnel region and 
stau-coannihilation are also shown to exist. With fiM2 > 0, the SUSY contributions 
to the muon anomalous magnetic moment can help provide better agreement than 
the SM with the experimental data. Finally, in comparison to SO(IO) with /i > 
and universal gaugino masses, there exist some important differences, even though 
the number of fundamental parameters in the two cases are the same. The lack of 
WMAP compatible neutralino dark matter in the yU > case is one of them. Also, 
with yU < 0, we find examples in which the first two squark families are relatively 
light (~ TeV), and the third family b and t squarks can be lighter than the gluino 
(which happens to be the lightest colored sparticle in SUSY SO(IO) with /x > ). 

Note Added: As we were finishing this work, Stuart Raby pointed out that 
M. Badziak, M Olechowski and S. Pokorski are also investigating SO(IO) Yukawa 
unification with /i < 
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